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ABSTRACT 


The  molecular  field  coefficients  employed  in  the  Neel  theory  of  ferrimagnetism 
have  been  determined  as  functions  of  the  levels  of  diamagnetic  ion  substitution 
in  the  garnet  family  {YzGd^_z)  [R^e^^]  (Qy.Fe.j_y)  ^12  w^ere  ^  an<*  Q  rep¬ 
resent  diamagnetic  octahedral  and  tetrahedral  substitutions,  respectively.  The 
coefficients  may  be  listed  as 


Ndd= 

-30.4 

a 

—  0.43x) 

N  = 
aa 

1 

O' 

© 

a 

-0.42y) 

N  = 

0 

cc 

Nad  = 

97.0  (1  - 

0.125x  —  0.127y) 

Ncd  = 

6.0 

Nac  = 

-3.44 

moles/cm 

With  these  coefficients  the  magnetic  moment  versus  temperature  curves  of 
compositions  ranging  from  0  <  x  <  0.70,  0^y^l.95,  and  0.40  ^  z  <  1 .00  were 
computed  and  are  presented  in  this  report. 
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MAGNETIC  MOMENT  VERSUS  TEMPERATURE  CURVES 
OF  FERRIMAGNETIC  GARNET  MATERIALS 


I.  INTRODUCTION 

The  phenomenon  of  ferrimagnetism  has  been  the  subject  of  extensive  investigation  since  its 
discovery.  From  a  phenomenological  standpoint,  the  molecular  field  theory  of  Neel  has  proven 

A 

to  be  extremely  successful  in  explaining  many  of  the  basic  properties  of  ferrimagnetic  systems. 

The  essence  of  the  simplest  model  consists  of  two  alternating  sublattices  of  unequal  and  anti¬ 
parallel  moments,  with  three  molecular  field  coefficients  employed  to  describe  the  exchange 
field  effects:  one  ferromagnetic  coefficient  for  each  sublattice  and  a  third  for  the  antiferromag¬ 
netic  interaction  between  the  sublattices.  Through  proper  use  of  Brillouin  functions  together 
with  the  correct  values  of  the  coefficients,  it  has  been  demonstrated  that  magnetization-temperature 
curves  may  be  accurately  reproduced.  Conversely,  experimental  curves  of  this  type  may  be 

used  to  determine  the  molecular  field  coefficients,  as  was  demonstrated  by  Anderson  using  a 

2 

trial-and- error  fit  for  yttrium-iron  garnet. 

In  many  cases,  the  coefficients  have  been  determined  from  magnetic  susceptibility  measure - 

3 

ments  above  the  Curie  temperature.  Measurements  of  this  type  were  reported  by  Pauthenet 

and  Aleonard  for  a  variety  of  iron  garnets.  A  third  method  involving  an  analytical  technique 

was  demonstrated  by  Rado  and  Folen  for  magnesium  and  lithuim  ferrite.^  In  each  of  these  cases, 

the  material  may  be  considered  as  an  ideal  Neel  ferrimagnet  (i.e.,  no  sublattice  canting)  since 

the  magnetic  moments  at  cryogenic  temperatures  approach  their  theoretical  values. 

3+ 

When  Fe  ions  are  replaced  by  diamagnetic  substitutes,  there  is  considerable  experimental 
evidence  that  sublattice  canting  or  some  equivalent  effect  begins  to  take  place.  In  a  model  de¬ 
veloped  by  Gilleo/*  departures  from  the  Neel  model  were  explained  by  considering  that  ions  be¬ 
came  paramagnetic  centers  when  interacting  with  less  than  three  nearest-neighbor  cations  of 
different  coordination.  For  small  substitutions,  this  model  proved  reasonably  satisfactory  but 
broke  down  quickly  at  higher  levels. 

7 

In  subsequent  work,  Geller  and  his  colleagues  performed  extensive  magnetization  measure¬ 
ments  on  substituted  yttrium -iron  garnets  at  cryogenic  temperatures  and  concluded  that  the 

g 

canting  originally  proposed  by  Yafet  and  Kittel  does  occur  and  in  a  peculiar  manner.  When 
substitutions  are  made  into  one  sublattice,  random  canting  takes  place  in  the  opposite  one.  As 
the  substitution  level  is  increased,  the  opposite  sublattice  gradually  reaches  an  antiferromagnetic 
state.  These  results  led  to  the  concept  of  competing  ferrimagnetic  and  antiferromagnetic  solvents. 
For  small  substitutions,  the  ferrimagnetic  solvent  dominates  the  antiferromagnetic  solvent;  after 
a  threshold  is  attained  at  some  higher  level  of  substitution,  the  reverse  becomes  true,  with  the 
result  that  an  abrupt  change  in  the  degree  of  canting  is  observed. 


1 


In  the  work  to  be  reported  here,  the  molecular  field  coefficients  of  several  yttrium-iron 

garnet  compositions  were  determined  from  the  experimental  magnetization-temperature  curves 
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reported  by  Gilleo  and  Geller.  The  results  reveal  unexpected  linear  relationships  between  the 
coefficients  and  levels  of  substitution  and  strongly  suggest  that  canting  is  directly  related  to  the 
changes  in  the  field  constants.  As  a  result,  an  attempt  is  made  to  establish  a  transition  between 
the  ideal  ferrimagnet  of  Neel  and  the  canted  sublattice  version  discussed  by  Geller,  as  the  amount 
of  diamagnetic  substitution  is  increased  from  zero  to  the  critical  regions  where  antiferromag¬ 
netism  becomes  dominant.  Within  the  ferrimagnetic  region  of  compositions,  it  should  now  be 
possible  to  compute  with  good  accuracy  the  magnetic  moment-temperature  curves  of  yttrium- 
iron  garnet  with  combinations  of  both  octahedral  and  tetrahedral  substitutions. 


H.  THEORETICAL  CONSIDERATIONS 
A.  N6el  Theory  of  Ferri magnetism 

According  to  the  Neel  model  of  ferrimagnetism,  the  temperature  dependence  of  the  magnetic 
moment  per  mole  of  each  sublattice  may  be  represented  by  Brillouin  functions 

Mt(T)  =  Mt(0)  Bg  (x)  ,  (1) 

i 

where  the  subscript  i  refers  to  the  particular  sublattice.  For  the  general  case  of  the  garnet 

structure,  three  sublattices  must  be  considered:  the  tetrahedral  (d),  the  octahedral  (a),  and  the 

dodecahedral  (c),  with  the  last  containing  rare-earth  ions.  In  this  work,  the  discussion  will  be 

3+ 

confined  to  the  two- sublattice  case,  i.e.,  d-  and  a-sites  containing  Fe  magnetic  ions.  The 
three -sublattice  case  is  outlined  in  Appendix  A. 

For  this  case,  the  magnetic  moment  per  mole  is  given  by 

M(T)  =  Md(T)  -  Ma(T)  ,  (2) 


where 

Md(T)  =  Md(0)  B  (xd) 
d 

Ma(T)  =  Ma(0)  Bs  (xa)  . 

a 

The  Brillouin  functions  are  expressed  as 
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In  Eq.  (4),  Ndd,  N&a,  and  Nda  =  Nad  are  the  molecular  field  coefficients,  and  S&  are  the  spin 
quantum  numbers  (5/2  for  the  high-spin  state)  of  the  Fe3+  ions  occupying  the  d-  and  a-sites, 
respectively,  g  is  the  spectroscopic  splitting  factor  (=  2.0),  is  the  Bohr  magneton,  and  k  is 
the  Boltzmann  constant. 

At  T  =  0°K,  the  magnetic  moments  per  mole  for  each  sublattice  are  given  by 

Md(0)  =  3gSd{xBN(l  -kd) 

Ma(0)  =  2gsaHLBN(1  "V  '  (5) 

3  + 

where  kd  and  k&  represent  the  fractions  of  diamagnetic  ions  substituted  for  Fe  ions  in  the  re¬ 
spective  sublattice  and  N  is  Avogadro's  number.  The  factors  3  and  2  appearing  in  Eq.  (5)  rep¬ 
resent  the  relative  numbers  of  d-  and  a-sites  in  the  garnet  formula  unit.  The  above  theory 
represents  the  essence  of  the  Neel  model,  which  is  based  on  antiparallelism  between  the  moments 
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of  the  two  sublattices.  As  indicated  in  Figs.  1  and  2,  which  are  replots  of  Geller' s  data,  there 
are  discrepancies  between  the  experimental  results  and  the  Neel  model.  For  kd  ^  0.65  and 
ka<  0.35,  the  differences  are  minor.  However,  above  these  values,  the  departures  become 
enormous  and  may  be  explained  only  by  assuming  significant  changes  in  the  spin  arrangments 
of  at  least  one  of  the  sublattices. 

According  to  Geller,  these  departures  from  the  Neel  model  are  caused  by  random  canting 
within  the  individual  sublattices.  Initially,  the  canting  is  relatively  small,  but  increases  quickly 
as  the  transition  points  are  reached  and  the  antiferromagnetic  tendencies  dominate.  In  the  work 
to  be  described,  it  will  be  necessary  to  compute  Md(0)  and  Ma(0)  below  the  transition  points,  and 
it  is  clear  from  Figs.  1  and  2  that  the  Neel  model  given  by  Eq.  (5)  will  not  be  adequate.  As  men¬ 
tioned  earlier,  Gilleo^  has  attempted  to  solve  this  problem  by  means  of  a  statistical  model.  Al¬ 
though  his  results  provided  better  agreement  with  the  data,  the  fit  was  not  satisfactory  over  the 
entire  range  of  interest  in  this  work,  i.e.,  k^  0.65  and  ka  ^  0.35.  In  order  to  reproduce  from 
theory  the  experimental  values  of  the  sublattice  moments  at  T  =  0°K  the  following  modifications 
to  Eq.  (5)  were  determined  empirically  as  the  only  convenient  recourse: 

Md(0)  =  3^VbN(1  “  kd}  (1  -  °*1  ka) 

Ma(0)  =  2gSaHBN(l  -  ka)  (1  -  k*’4)  (6) 


There  is  no  obvious  explanation  of  why  a  linear  factor  appears  to  work  in  one  case  while  an 
exponential  will  suffice  in  the  other.  Clearly,  both  factors  are  but  approximations  to  a  more 
general  theory  yet  to  be  devised.  As  will  be  demonstrated  later,  the  answer  is  probably  related 
to  the  initial  magnitudes  of  the  intra- sublattice  molecular  field  coefficients  or  the  relative  strengths 
of  the  d  and  a  sublattice  exchange  energies. 

Before  this  section  is  completed,  it  should  be  mentioned  that  all  of  the  magnetic  moments 
in  the  figures  to  follow  are  expressed  in  Bohr  magnetons  per  formula  unit,  or 


V 


M 

B'  ^bn 


(7) 


where  pertinent.  All  other  quantities  involved  were  expressed  in  units  of  the  cgs  system. 
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B.  Theoretical  Model 

The  basic  approach  to  this  problem  was  the  trial -and -error  determination  of  the  correct 
combination  of  N^,  N&a,  and  required  to  fit  a  given  set  of  data.  For  pure  Y^Fe^O^,  a 
set  of  coefficients  had  already  been  furnished  by  Anderson.^  Since  the  molecular  field  coefficient 
is  directly  related  to  the  strength  of  an  exchange  field,  it  is  not  possible  to  justify  increases  in 
magnitude  with  diamagnetic  substitutions.  As  a  result,  only  reductions  in  the  magnitudes  of  the 
coefficients  were  considered  as  physically  realistic.  It  also  became  apparent  that  when  substi¬ 
tutions  were  confined  to  one  sublattice,  it  was  not  possible  to  obtain  a  good  fit  with  experiment 
by  reducing  the  magnitudes  of  both  and  Naa  simultaneously,  because  excessive  reductions 
in  the  Curie  temperature  led  to  increases  in  From  these  observations,  the  following  con¬ 

clusions  were  drawn:  (1)  to  a  first  approximation,  substitutions  in  the  d  sublattice  cause  reduc¬ 
tions  in  | Naa |  without  changing  |N^|,  and  vice  versa,  and  (2)  |Na^|  decreases  regardless  of 
the  site  distribution  of  the  diamagnetic  ions. 

To  place  this  reasoning  on  a  firmer  physical  basis,  consider  ar\  ion  of  moment  fT  in  a  molecu¬ 
lar  field  IT  (Ref.  10).  The  interaction  energy  may  be  expressed  as 


E=~tV  H. 


(8) 


E  —  VBV  2  V,  ' 


(9) 
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Since  this  energy  is  equivalent  to  the  sum  of  the  exchange  interactions  with  its  near  neighbors, 

(10) 


E  =  -2  y  Z..J..S.  •  S.  , 
u  13  13  i  3 

3 


.th 


where  z..  represents  the  number  of  nearest  neighbors  on  the  j  sublattice  that  interact  with  the 
th 

i  ion  and  is  the  exchange  constant.  Since 


M.  =  n.g.iiT-jS.  , 
3  3*3^B  3 


(11) 


where  n.  is  the  number  of  ions  per  mole  in  j*h  sublattice,  the  molecular  field  constant  may  be 
found  by  combining  Eqs.  (9),  (10),  and  (11): 


J.. 


N..  =  (-ii)2  _ii, 
1J  '  "j  '  gjg^B 


(12) 


In  cases  where  random  substitutions  are  considered,  (z„/m)  is  unaffected  (assuming  z_  is  treated 
as  an  average)  and  the  only  variable  is  J^,  which  must  also  be  considered  as  an  average.  There¬ 
fore,  the  molecular  field  coefficients  should  be  directly  proportional  to  the  average  exchange 
constants  as 

Ndd a  <Jdd> 

Naa  “  <Jaa> 


Nad  a  <Jad> 


(13) 
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When  a  magnetic  ion  is  removed  from  the  d  sublattice,  as  sketched  in  Fig.  3,  all  of  the 

nearest -neighbor  a  cations  experience  the  loss  of  a  fraction  of  their  superexchange  interactions 

through  canting.  In  other  words,  the  local  J  and  J  are  weakened  to  J'  ,  J"  ,  and  J’  ,,  while 

aa  a  a  aa  aa  ad 

to  a  first  approximation,  the  d  sublattice  is  affected  only  insofar  as  it  suffers  the  loss  of  moment 
associated  with  the  missing  ion. 

Both  sublattices  consequently  suffer  reductions  in  molecular  field  strengths,  one  through 
reduced  exchange  and  the  other  through  a  loss  in  total  magnetic  moment.  However,  as  indicated 
by  Eq.  (13),  only  N  and  N  ,  will  change  because  of  lower  values  of  <J  >  and  <J  Since 

3.3  aU  3.3.  30 

<Jdd^  *s  unaffected,  remains  the  same.  These  general  arguments  also  apply  to  the  case 
where  an  a-site  magnetic  ion  is  removed,  with  the  result  that  N&a  will  be  unaffected. 

The  relationship  of  this  model  to  that  of  Geller  is  almost  axiomatic.  If  the  molecular  fields 
are  considered  as  bases  for  the  parallelism  of  the  sublattices,  it  is  evident  that  any  tendency 
toward  canting  would  of  necessity  depend  on  reductions  in  the  magnitudes  of  these  coefficients. 
Since  Geller  has  pointed  out  the  tendency  toward  antiferromagnetism  at  higher  substitutions,  it 
seems  likely  that  this  trend  should  initially  be  manifested  by  changes  in  the  values  of  N^,  N&a, 
and  This  subject  will  be  discussed  more  fully  in  a  later  section. 

IE.  DETERMINATION  OF  MOLECULAR  FIELD  COEFFICIENTS 

In  Figs.  4,  5,  8,  and  9,  the  ’’best-fit"  magnetic  moment  versus  temperature  curves  of  several 

q 

compositions  are  plotted  together  with  the  experimental  points  determined  by  Gilleo  and  Geller. 
For  purposes  of  comparison,  the  dashed-line  curves  computed  by  using  the  molecular  field  co¬ 
efficients  of  pure  Y-jFe^O^  (see  Table  I)  are  also  included.  The  sequence  employed  in  this 


TABLE  1 

MOLECULAR  FIELD  COEFFICIENTS  OF  SUBSTITUTED  Y0FecO 

0  0 

12 

Ndd 

N 

oa 

Z 

a 

Q_ 

Composition 

kd 

k 

a 

moles/cm^ 

moles/cm^ 

moles/cm 

Y3Fe5°12 

0 

0 

-30.4 

-65.0 

97.0 

Y3Sc0.  25Fe4.  75°  12 

0 

0.125 

-27.0 

-65.0 

94.0 

Y3ln0.50Fe4. 50°12 

0 

0.25 

-24.0 

-65.0 

91.0 

Y3Go0.25Fe4.75°12 

0.075 

0.0125 

-30.1 

-59.0 

94.0 

Y3Ga0.75Fe4.25°12 

0.22 

0.04 

-29.4 

-47.0 

88.0  , 

Y3Al0.33Fe4.67°12 

0.105 

0.007 

-30.2 

-56.0 

93.0 

Y3Al1.0Fe4.0°12 

0.28 

0.07 

-28.6 

-42.0 

85.0 
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3+  3  + 

work  was  as  follows:  from  pure  octahedral  (a-site)  substitutions  of  Sc  and  In  .shown  in  Figs.  4 
and  5,  the  curves  were  fitted  by  reducing  |Ndd|  and  |Nad|.  The  results  were  then  plotted  as 
functions  of  ka  as  in  Figs.  6  and  7,  where  it  was  found  that  linear  relations  exist  in  both  cases. 
From  the  slopes  of  these  straight-line  curves, 

N  ,  ,  =  -30.4  +  26.4k 
dd  a 


Nad  =  97,0  ”24-0ka  (for  kd  =  0)  *  (14) 

3+  3  + 

With  these  expressions  it  was  then  possible  to  fit  the  results  for  the  Ga  and  Al  substituted 
compositions  shown  in  Figs.  8  and  9.  After  determining  the  values  of  k&  and  kd  from  the  moments 
at  T  =  0°K,  the  values  of  Ndd  for  all  four  of  these  compounds  were  calculated  from  Eq.  (14). 

These  values  were  then  used  in  determining  the  proper  theoretical  fits,  which  were  obtained  this 
time  by  varying  Naa  and  Nad.  The  results  for  the  molecular  field  coefficients  of  all  compositions 
studied  are  listed  in  Table  I. 

As  before,  the  values  of  Naa  were  found  to  have  a  linear  dependence  on  k^,  as  indicated  by 
the  appropriate  curve  in  Fig.  6.  In  the  case  of  N  the  results  from  Fig.  7  were  employed  with 
(Na(j  +  24  ka)  plotted  as  a  function  of  k^  in  Fig.  10  in  order  to  reveal  the  linearity  and  determine 
the  correct  slope.  The  results  may  thus  be  expressed  as 

N  ,  ,  =  -30.4  (1  -0.87k  ) 
dd  a 

N  =  -65.0  (1  -1.26k.) 
aa  d 

Nad  =  97,0  (1  — °.25ka“°-38kd)  '  (15) 

subject  to  the  limitations  of  k  ^0.35  and  k ,  ^  0.65. 

a  d 

Since  kd  =  y/3  and  ka  =  x/2,  where  y  and  x  are  the  numbers  of  substituted  diamagnetic  ions 
per  formula  unit  occupying  d-  and  a-sites  respectively,  Eq.  (15)  may  be  written  as 

Ndd  =  -30.4  (1  —  0.43x) 


Naa  =  -65.0  (1  -  0.42y) 

Nad  =  97,0  (1  “  °*125x  “  °-127y)  (16) 


for  x<  0.70  and  y^  1.95 

AN  , 
ad 


N;d 


From  Eq.  (16)  it  may  be  readily  deduced  that 
/IAN..  AN  |\ 

+0,3(hd +  bin)  * 


(17) 


where  N*  =  -30.4,  N*  =-65.0,  and  N°.  =  97.0. 
dd  aa  ad 

It  is  indeed  curious  that  the  numerical  coefficients  of  x  and  y  in  Eq.  (16)  are  so  similar. 

Although  there  is  no  explanation  to  be  offered  at  this  time,  the  following  observation  will  be 

3+  3  + 

made.  When  a  d-site  Fe  ion  is  replaced,  four  nearest-neighbor  a-site  Fe  ions  are  directly 

affected.  However,  when  the  reverse  occurs,  six  nearest  neighbors  are  involved,  an  increase 

of  fifty  percent.  Thus,  the  ratio  of  the  number  of  sites  affected  to  the  total  number  of  sites 

within  a  given  sublattice  is  effectively  the  same  in  both  cases.  Therefore,  it  must  follow  that 

the  contributions  to  the  percentage  decrease  in  the  coefficient  from  each  affected  site  must  also 

be  the  same  in  both  cases,  in  order  to  produce  the  above  relations  for  Ndd  and  Naa«  As  a  direct 

result  of  this,  it  should  not  be  surprising  to  find  the  symmetry  expressed  in  Eq.  (17). 
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With  the  results  given  by  Eq.  (15),  it  should  now  be  possible  to  compute  the  magnetic  moment- 
temperature  curve  for  any  substitution  below  the  antiferromagnetic  transition  points.  In  Figs.  11 
and  12,  normalized  curves  are  plotted  for  pure  tetrahedral  and  octahedral  substitutions  respec¬ 
tively,  and  reveal  the  extent  to  which  the  shape  of  the  curve  may  vary  as  the  molecular  field 
coefficients  change.  As  a  test  for  systems  containing  diamagnetic  ions  in  both  sites,  the  ex¬ 
perimental  results  for  {Y^}  [MgxFe2xl  (Fe3_xSi  )  012  (Ref.  7)  are  compared  with  theory  in 
Fig.  13.  The  agreement  is  excellent  up  to  the  levels  of  substitution  indicated,  without  any  indi¬ 
cation  of  a  possible  breakdown  at  higher  concentrations. 

IV.  SUBLATTICE  CANTING 

In  the  development  of  Eq.  (6),  it  was  stated  that  random  canting  of  the  sublattice  opposite  to 
that  in  which  substitutions  are  made  was  concluded  from  the  extensive  experimental  work  of 

7 

Geller.  Since  canting  effects  were  observed  more  readily  in  the  tetrahedral  sublattice,  it  was 
decided  that  its  tendency  towards  antiferromagnetic  alignment  was  greater  than  for  the  octahedral 
case.  In  this  sense,  d-d  interactions  were  considered  to  be  stronger  than  a-a  interactions.  In 
this  present  work,  the  pseudo-ferromagnetism  of  the  individual  sublattices  is  chosen  as  the  basis 
for  discussion  because  of  its  obvious  relation  to  the  sign  and  magnitudes  of  N^  and  Na&.  From 
this  standpoint,  the  tetrahedral  sublattice  is  considered  to  have  the  weaker  ferromagnetic  inter¬ 
action  because  In,, I  <  |N  |. 

1  dd  aa 

In  the  model  depicted  in  Fig.  3,  the  nearest-neighbor  magnetic  moments  of  the  sublattice 
are  shown  as  canted  when  a  particular  d  ion  is  removed.  If  indeed  this  occurs,  it  would  be 
surprising  to  find  that  a  reduction  in  |N  |  did  not  also  occur.  In  Fig.  14,  the  percentage  change 
from  canting  in  the  Bohr  magnetons  per  formula  unit  at  T  =  0°K,  £17^/17^,  is  plotted  as  a  function 
of  the  sublattice  molecular-field  coefficient  N„.  The  results  point  to  a  close  relationship  be¬ 
tween  these  two  variables,  as  both  curves  approach  the  same  linear  function  just  prior  to  the 

3 

antiferromagnetic  transition  which  sets  in  at  N^  =*  —  20  moles/cm  , 

At,B  -3 

— p  = -(3.75NU  +  11.4)  X  10  .  (18) 

nB 

As  suggested  by  Geller,  the  canting  effects  theoretically  should  be  expected  immediately 
upon  diamagnetic  substitutions  into  either  sublattice.  For  the  d  sublattice,  the  change  is  im¬ 
mediately  detectable,  as  seen  by  Fig.  2,  because  of  its  relatively  weak  molecular  field  coefficient 
N^  =  —30.4.  With  Naa  =  —65.0,  a  significant  amount  of  tetrahedral  substitution  is  required  be¬ 
fore  the  coefficient  drops  to  the  point  where  departures  from  the  Neel  model  are  detected. 

V.  CONCLUSIONS 

In  the  foregoing  sections,  a  description  of  the  method  used  in  determining  the  molecular 

field  coefficients  of  several  substituted  yttrium-iron  garnets  was  presented.  The  coefficients 

were  found  by  fitting  the  Neel  theory  to  experimental  magnetic  moment  versus  temperature 

curves  below  the  antiferromagnetic  transition  substitution  levels  of  k^<  0.65  and  ka  <  0.35. 

In  general,  it  was  found  that  peculiar  linear  relations  exist  between  the  coefficients  N^, 

N  ,  and  N  ,  and  the  fractions  of  sublattice  substitution  k  ,  and  k  .  For  substitutions  in  one 
aa  ad  da 

sublattice,  the  coefficient  of  the  opposite  sublattice  is  reduced  in  magnitude  while  its  own  co¬ 
efficient  remains  unchanged;  for  substitutions  in  either  sublattice,  the  inter -sublattice  coefficient 
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is  reduced.  By  comparing  the  magnitudes  of  the  coefficients  with  the  changes  in  sublattice 
moment  from  canting,  it  was  found  that  at  these  lower  levels  of  substitution,  the  amount  of  canting 
in  a  given  sublattice  appears  to  be  directly  related  to  a  change  in  its  molecular  field  coefficient. 
This  result  provides  a  tie  between  the  Neel  theory  and  the  random  canting  concepts  of  Geller 
and  serves  to  confirm  that  canting  begins  to  set  in  immediately  upon  substitution. 

Since  |Ndd|  <  |NagJ  initially,  canting  would  be  expected  to  occur  earlier  in  the  tetrahedral 
sublattice,  in  accord  with  Geller^  observations.  This  result  also  conforms  with  the  notion  that 
the  d  sublattice  has  a  stronger  antiferromagnetic  interaction,  because  the  weaker  |N^^|  co¬ 
efficient  refers  to  the  strength  of  ferromagnetic  interaction.  If  it  is  weaker  ferromagnetically, 
it  must  follow  that  it  is  stronger  antiferromagnetically. 

Finally,  because  it  is  possible  to  express  the  coefficients  as  analytical  functions  of  the 
amounts  of  substitutions,  it  appears  likely  that  for  ions  of  the  valence  states  involved  in  this 
study,  magnetic  moment-temperature  curves  may  be  computed  readily  for  any  level  of  substi¬ 
tution  below  the  antiferromagnetic  transitions.11  In  situations  where  the  levels  of  substitution 
in  both  sublattices  approach  their  limits,  some  question  of  applicability  still  remains.  However, 
the  results  shown  in  Fig.  13  suggest  that  there  is  good  reason  for  optimism  on  this  point. 
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APPENDIX  A 


For  the  three  sublattice  case,  the  Neel  theory  must  be  augmented  as  follows: 


where 


M  =  M  +  M  -M  ,  , 

a  c  d  ' 


M.(T)  =  M.(0)  B  (x.) 

1  1  Oj  1 

s.gnR 

x,  =  (N.  ,M  ,  +  N,  M  +  N,  M  ) 

d  kT  dd  d  da  a  dc  c 

x  =  -4^  (N  ,M ,  +  N  M  +  N  M  ) 
a  kT  ad  d  aa  a  ac  c 

S  gp.-p. 

x  =  -g-yfr-  (N  ,M  -  +  N  M  +  N  M  ) 
c  kT  cd  d  ca  a  cc  c 

with  N  ,  =  N  ,,  N,  =  N  ,  and  N  =  N  . 
da  ad*  dc  cd  ac  ca 

At  T  =  0°K,  the  magnetic  moments  per  formula  unit  are  given  by 
Md(0)  =  3gSdPBN(l  -kd)  (1  -0.1ka) 

Mc(0)  =  3gSctxBN(l  -  kc)  . 


(A-l) 


(A-2) 


(A-3) 


3+  o  _j_ 

where  k  is  the  fraction  of  Y  ions  in  c -sites  and  S  =  7/2  for  Gd  . 

^  3+  ^  ^ 

For  Gd  substitutions  into  the  c- sublattice,  the  Neel  model  is  used  without  modification 

because  the  effects  of  canting  both  in  the  sublattice  itself  and  in  the  other  sublattices  are  judged 

to  be  negligible  at  the  levels  of  substitution  involved.  This  conclusion  is  based  on  inspection  of 
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the  work  of  Geller  et  al.,  who  reported  data  on  the  garnet  families  containing  Gd  ions.  To 

a  first  approximation,  the  molecular  field  coefficients  of  the  c- sublattice  are  considered  to  be 

independent  of  composition. 

The  molecular  field  coefficients  may  be  listed  as  follows: 

Ndd  =  -30.4  (1  -  0.87  ka) 

N  =  -65.0  (1  -  1.26k.) 
aa  a 

=  0 
cc 

In  Figs.  15,  16,  and  17,  comparisons  between  theory  and  experiment  are  presented  to  illus¬ 
trate  the  accuracy  of  the  computed  curves.  Those  discrepancies  which  may  not  be  accounted 
for  by  experimental  error  are  probably  the  result  of  the  small  canting  effects  which  have  been 
ignored.  Since  the  data  in  Fig.  17  were  not  corrected  for  porosity  (which  could  have  been  as 
high  as  5  percent),  it  may  be  assumed  that  the  agreement  could  be  significantly  improved. 


N  .  =  97.0  (1  -  0.25k  -0.38k.) 
ad  a  d 


N  ,  =  6.0 
cd 


N  =  —  3.44moles/cm“ 
ac  ' 


(A-4) 
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APPENDIX  B 
COMPUTER  PROGRAM 
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C  MAIN  PROGRAM  FOR  3  MOMENTS--  USING  SUBROUTINE  M0M3  FOR  CON- 
C  VERGENCE.  G.  DIONNE 

IMPLICIT  REAL-8  (A-H,K,M,N,0-Z) 

COMMON  /Ml/  SD, SA, SC,MUB, K, BN, G 
COMMON  /M2/  NDD,NAA,NAD,NAC,NCD,NCC 
COMMON  /M3/  KD, KA, KC 
COMMON  /ANS/  M/MD1/MA1/MC1/ NBANS 

DIMENSION  TEMP(200), NB( 200 ), KAA( 21 ), KDD( 21),NBB(21, 200 ) ,MSAVE( 21 ), 
1KCC( 5 ) 

C  THE  FOLLOWING  ARRAYS  ARE  USED  IN  PLOTTING  ONLY 
DIMENSION  TITLE(9)/YY(1700)/TT(1700)/LC(21) 

C  T1  AND  T2  ARE  LABELS  FOR  THE  X  AND  Y  AXES 
DATA  T1/,TEMP'/,T2/'NB  '/ 

DATA  KCC/.8,.7,.6,.5,.4/ 

THE  FOLLOWING  TWO  LINES  REFER  TO  ROUTINES  ACCESSIBLE  TO  USERS  OF 
LINCOLN'S  360/67  AND  PERTAIN  ONLY  TO  THE  PLOTTING 
CALL  REREAD( 15, 80) 

CALL  STOIDVC  G.  DIONNE',9,0) 

SD-2.5 
SA-2.5 
SC-3.5 

MUB-9. 27E-21 
K-1.38E-16 
BN-6. 023E+23 
G-2.0 
NCD-6.0 

NAC— .  5734-NCD 
NCC-0.0 

INITIAL  KA  VALUE 
KAA( l)-.05 
INITIAL  KD  VALUE 
KDD( 1)-. 35 

INCREMENT  FOR  KA  AND  KD 
DK-.01 

INCREMENT  FOR  TEMPERATURE 
DT-10. 

NUMBER  OF  CURVES  PER  PAGE  (NUMBER  OF  KD'S  USED) 

I  CURVE-11 
DO  59  L-l, I  CURVE 
KAA( L) -KAA( l)+( L-1)*DK 
59  KDD(L)»KDD(1)+(L-1)*DK 
DO  218  LRR-1,5 
KC-KCC(LRR) 

KC-1.0 

DO  100  LKA-1,6 
KA-KAA(LKA) 

C  THE  FOLLOWING  3  LINES  ARE  RELEVANT  ONLY  TO  THE  PLOTTING 
IXZ-0 

NMAX  — 25.E+25 
NMIN-25.E+25 
DO  200  LKD-1, (CURVE 
KD-KDD(LKD) 

NDD--30. 4*( 1. -. 87*KA) 

NAA--65. *(1.-1.26*KD) 

NAD-97. *(1.-.25*KA-.38*KD) 

TEMP(l)-0.0 

TEMP(2)-20. 

IFCKD.LT.. 10.0R.KA.LT.. 10. OR.KC.LT.. 10)  TEMP(2)-40. 

IF(TEMP(2).EQ,40.)  GO  TO  21 

TEMP(29)-293. 

TEMP(35)-343. 

DO  31  L-2,28 

31  TEMP(L)-TEMP(2)+(L-2)*DT 
DO  32  L-30,34 
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32  TEMP(L)-TEMP(28)+(L-29)*DT 
DO  33  L-36,72 

33  TEMP(L)-TEMP(34)«-(L-35)*DT 
GO  TO  22 

21  TEMP(28)-293. 

TEMP(34)«343. 

DO  34  L«2,27 

34  TEMP(L)-TEMP(2)*(L-2)*DT 
DO  35  L-29, 33 

35  TEMP(L)-TEMP(27)*(L-28)*DT 
DO  36  L»35,72 

36  TEMP(L)-TEMP(33)+(L-34)*DT 

22  CALL  M0M3 
Ml-1. 

M2— 9999. 

DO  70  L-l, 72 
T-TEMP(L) 

CALL  C0NTIN(T, ITER, IER) 

IF(IER.EQ.l)  GO  TO  150 
IF(TEMP(2).EQ.20. .AND. L. EQ. 29 )  Ml-M 
IF(TEMP(2).EQ.20..AND.L.EQ.35)  M2-M 
IF(TEMP(2).EQ.40. .AND. L.EQ.28)  Ml-M 
I F (TEMP( 2 ) . EQ. 40 . .AND. L. EQ. 34 )  M2-M 
NBB(LKD,L)«NBANS 
70  CONTINUE 

150  IF(IER.EQ.l)  LC(LKD)-L-1 
IF(IER.EQ.O)  LC(LKD)-L 
MSAVE( LKD) -(M1-M2 )/Ml 
LE-LC(LKD) 

THE  FOLLOWING  DO  LOOP  FORMS  THE  ARRAY  OF  ABSCISSAS  FOR  THE  POINTS 
TO  BE  PLOTTED 
DO  77  LA-1, LE 
IXZ-IXZ+1 

77  TT( IXZ)-TEMP(LA) 

200  CONTINUE 

WRITE(6,4)  KC 
4  FORMAT(//'  KC-',F5.2) 

WR I TE( 6, 888 ) 

DO  98  LN-1, I  CURVE 

98  WRI TE(6,999)  KA, KDD( LN), NfiB( LN, 1 ), TEMP( LC( LN) ),MSAVE( LN) 

888  FORMAT (/10X, ' KA' , 7X, '  KD ' , 7X, '  NB  AT  TEMP-0 ', 7X, ' TEMP  LIMIT', 

17X, 1 (M1-M2)/M1' ) 

999  FORMAT (8X,F5.3,4X,F5.3,8X,F8.3,10X,F8.3,7X, E12.5) 

C  THE  FOLLOWING  DO  LOOP  FORMS  THE  ARRAY  OF  ORDINATES  FOR  THE  POINTS 
C  TO  BE  PLOTTED 
IN-0 

DO  61  LKD-1, I  CURVE 
LP-LC(LKD) 

DO  61  LZ-1, LP 
IN-IN+1 

I F ( LC( LKD) . EQ. 1 )  GO  TO  61 
I F ( NBB( LKD, LZ) .GT. NMAX )  NMAX-NBB( LKD, LZ) 

I F( NBB( LKD, LZ) . LT. NMI N)  NMI N-NBB( LKD, LZ) 

61  YY( I N)-NBB( LKD, LZ ) 

C  THE  NEXT  6  LINES  WRITE  ON  A  9-TRACK  TAPE  THE  INFORMATION  WHICH 
C  DATAGRAPHIX'S  D-4060  CONVERTS  INTO  FINISHED  PLOTS 
WR I TE( 15,3)  KC, KA, ( KDD( LSD) , LSD-1,  I  CURVE) 

3  FORMAT ( '  KC- ' , F5. 3, '  KA»',F5.3,'  KD- ' , 10 ( F3 . 2, ' , ' ), F3. 2 ) 

READ( 15,39)  TITLE 
39  FORMAT ( 9 A8) 

CALL  GRAPH (TT,YY, I  CURVE, LC, TITLE, 2) 

CALL  LINEAR(0.,740.,0.,NMAX,T1,T2,2) 

100  CONTINUE 
218  CONTINUE 

C  THIS  CALL  SIGNALS  THE  END  OF  PLOTTING 
CALL  PLTND 

99  CALL  EXIT 
END 
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CALCULATION  OF  M  VS  TEMP  FOR  GARNETS  WITH  GD  (+03) 

SUBROUTINE  M0M3 

IMPLICIT  REAL*8  (A-H,K,M,N,0-Z) 

COMMON  /Ml/  SD,SA,SC,MUB,K,BN,G 
COMMON  /M2/  NDD, NAA, NAD, NAC, NCD, NCC 
COMMON  /M3/  KD, KA, KC 
COMMON  /ANS/  M,MD1,MA1,MC1, NB 
P-MUB*BN 

MDO-3.*G*SD*P*(l.-KD)*(l.-.l*KA) 

MAO-2. *G*SA*P*(1.-KA)*(1.-KD**(5.4)) 

MCO-3.*G*SC*P*(l.-KC) 

C1-1./(2.*SD) 

C2-(2.*SD«-1.)*C1 
C3-1./(2.*SA) 

C4-(2.*SA+1.)*C3 
C5-1./(2.*SC) 

C6-(2.*SC+1.)*C5 
MD-MDO 
MA-MAO 
MC-MCO 
MD1-MDO 
MAI “MAO 
MC1-MCO 
RETURN 

C*********************************************************************** 
ENTRY  CONTI N(T, I  TER, I ER) 

ITER-0 

IER-0 

IF(T.EQ.O.O)  M“DABS ( MD-MA-MC ) 

IF(T.EQ.O.O)  GO  TO  55 
FR-G*MUB/(K*T) 

IX-1 

9  ITER-ITER+1 

I F( ITER. GE. 300)  GO  TO  999 
XD-SD*FR*(NDD*MD«-NAD*MA+NCD*MC) 

XA-SA*FR*(NAA*MA«-NAD*MD«-NAC*MC) 

XC-SC*FR*(NCC*MC+NCD*MD+NAC*MA) 

BSD-C2*DCOSH(XD*C2 )/DS I NH(XD*C2 )-Cl*DCOSH(XD*Cl)/DS I NH(XD*C1) 

BSA»C4*DCOSH(XA*C4)/DSINH(XA*C4)-C3*DCOSH(XA*C3)/DSINH(XA*C3) 

BSC-C6*DCOSH(XC*C6)/DSI NH(XC*C6)-C5*DCOSH(XC*C5)/DSI NH(XC*C5) 

MD2“MDO*BSD 

MA2“MAO*BSA 

MC2-MCO*BSC 

GO  TO  (13,14), IX 

13  I X“2 
MD1-MD 
MA1-MA2 
MC1-MC2 
GO  TO  15 

14  IX-1 
MD1-MD2 
MA1-MA 
MC1-MC2 

15  I F (DABS(MAl-MA) . LT. . 0001. AND. DABS (MD1-MD) . LT..0001 
1. AND. DABS(MC1-MC).LT.. 0001)  GO  TO  8 

MD-MD1 

MA-MA1 

MC-MC1 

IFCMD.LT. 0.0. OR.MA.LT. 0.0. OR.MC.LT. 0.0)  GO  TO  999 
GO  TO  9 

8  M“DABS (MD-MA-MC) 

55  NB»M/(MUB*BN) 

GO  TO  888 
999  IER-1 
888  RETURN 
END 
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APPENDIX  C 


By  systematically  varying  parameters  kc,  and  ka,  magnetic  moment  versus  temperature 
curves  for  a  wide  range  of  diamagnetic  substitutions  were  generated  from  the  computer  program 
stated  in  Appendix  B.  Within  the  limits  established  in  the  text,  k&  was  varied  in  steps  of  0.05 
between  0  and  0.35.  The  parameter  k^  was  varied  between  0  and  0.50,  also  in  steps  of  0.05, 

except  for  k  =1.0  where  it  was  extended  to  0.65.  Finally,  the  parameter  k  ,  which  defines  the 

c  3+  3+  c 

relative  amounts  of  Y  and  Gd  in  the  dodecahedral  sites,  was  varied  between  0.4  and  1.0  in 

steps  of  0.10.  Values  of  kc  below  0.4  were  not  included  because  they  were  considered  to  repre¬ 
sent  compositions  outside  of  the  current  range  of  practical  interest.  If  desired,  they  may  be 
readily  calculated  following  the  procedure  outlined  in  the  computer  program. 

The  results  of  these  computations  are  presented  in  Figs.  18  through  25. 
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Fig.  1.  Replot  of  Geller's  magnetic  moment  data  at  T  =  0°K  for  pure  tetrahedral 

substitutions,  indicating  a  departure  from  the  N6el  model.  Transition  to  the 

antiferromagnetic  state  occurs  at  k  ,  “  0.65. 

d 


Fig.  2.  Replot  of  Geller*s  magnetic  moment 
data  at  T  =  0°K  for  pure  octahedral  substitu¬ 
tions,  indicating  a  departure  from  the  Neel 
model.  Transition  to  the  antiferromagnetic 
state  occurs  at  k  “0. 35. 
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Fig.  3.  Two-dimensional  model  of  ferrimagnetic  sublattices,  suggesting  possible 
effects  on  near-neighbor  magnetic  moments  when  a  diamagnetic  substitution  is 
made  in  the  d  sublattice.  Lower  part  indicates  that  canting  of  local  octahedral 
site  moments  may  occur  when  a  magnetic  ion  is  removed  from  a  tetrahedral  site. 
The  argument  applies  equally  well  to  tetrahedral  site  canting  from  octahedral 
site  substitutions. 
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Fig.  4.  Comparison  of  theory  and  experiment 
for  magnetic  moment-temperature  curve  of 
Y3Sco,25^e4.  7501 2-  Dashed  line  represents 
the  result  calculated  by  using  the  molecular 
field  coefficients  of  pure  Y3Fe5<D]2* 


T  CK  ) 


Fig.  5.  Comparison  of  theory  and  experiment 
for  magnetic  moment-temperature  curve  of 
YsIhO.  5^e4. 50l2#  Dashed  line  represents 
the  result  calculated  by  using  the  molecular 
field  coefficients  of  pure  Y3Fe50]2* 
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Fig.  6.  Linear  variations  of  and  Naa 
with  kQ  and  kj,  respectively.  Upper  curve 

•s  for  Naa. 


Fig.  7.  Linear  variation  of  Nqcj  with  kQ  (k^  =  0). 
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Fig.  8.  Comparison  of  theory  and  experiment 
for  magnetic  moment-temperature  curves  of 
Y3Gao.  25Fe4. 75°  12  and  Y3CX10. 75Fe4. 250 12- 
Dashed  lines  represent  the  results  calculated  by 
using  the  molecular  field  coefficients  of  pure 
Y3Fe5<Di2. 


T  CK) 


200  400  600 


T  CK) 


Fig.  9.  Comparison  of  theory  and  experiment 
for  magnetic  moment-temperature  curves  of 

Y3Al0.33Fe4.67°12  and  Y3Al]Fe40i2* 

Dashed  lines  represent  the  results  calculated 
by  using  the  molecular  field  coefficients  of 
pure  Y3Fe50i2. 
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Fig.  11.  Theoretical  normalized  variations  of  magnetic  moment  with  temperature 
for  pure  tetrahedral  site  substitutions. 
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Fig.  12.  Theoretical  normalized  variations  of  magnetic  moment 
with  temperature  for  pure  octahedral  site  substitutions. 


Fig.  13.  Comparison  of  theory  with  experiment  for  three  compositions 
of  {Y3}  [Mgx^e2-x^  ^e3-x^'x)  fea*ure  substitutions  in 

both  sublattices. 
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Fig.  14.  Percent  change  in  sublattice  magnetic  moment  (attributed  to  canting)  as  a  function 
of  the  sublattice  molecular  field  coefficient.  At  Njj “  —  20moles/cm3,  antiferromagnetic 
transitions  occur  in  both  cases. 


Fig.  15.  Comparison  of  theory  and  experiment  for  magnetization-temperature 

curves  of  Y^Gd^Fe^^  for  z  =  0.6,  1.2,  1.8,  2.4,  and  3.0. 
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Fig.  16.  Comparison  of  theory  and  experiment  for  magnetization-temperature 
curves  of  Y^Gd^Fe^O^  for  *  =  0. 3,  0. 9,  1.5,  2.1,  and  2.  7. 


Fig.  17.  Comparison  of  theory  and  experiment  for  magnetization-temperature 
curves  ofY^  ^Gdj  ^Al^Fe^^O^  for  x  =  0/  0*  1  #  0*  2,  and  0.  5. 
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Fig.  18(a-h). 


Magnetic  moment  n^ 


vs  T  fork  =1.0,  0<k,^0.65,  k  (as  indicated), 
c  da 
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Fig.  18(a-h).  Continued. 
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Fig.  18(a~h).  Continued. 
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Fig,  18(a-h).  Continued. 
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T  CK) 

Fig.  19.  Magnetic  moment  |nJ  vs  T  for  k  =  k  =  0,  0<  k  <1.0. 

d  a  a  c 
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(a)  For  k  =  0. 
x  7  a 


(b)  Fork  =0.05. 
x  a 


Fig.  20(a-h).  Magnetic  moment  |n^|  vs  T  for  k^  =  0.90,  0^  k^  ^  0.50 ,  k^  (as  indicated). 
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(c)  Fork  =0.10. 
N  '  a 
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(d)  For  k  =0.15. 
a 


Fig.  20(a-h).  Continued. 
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Fig.  20(a-h).  Continued. 
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(h)  For  kQ  =  0.35. 
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Fig.  20(a-h).  Continued. 
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(a)  For  k  =  0. 
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(b)  Fork  =0.05. 
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Fig.  21(a-h).  Magnetic  moment  InJ  vs  T  fork  =  0. 80,  0  k  4  0. 50,  k  (as  indicated). 

b  e  da 
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Fig.  21(a-h).  Continued. 
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(e)  For  k  =  0. 20. 
'  a 
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Fig.  21(a-h).  Continued. 
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(g)  For  k  =  0. 30. 
9  a 


(h)  For  kQ  =  0. 35. 


Fig.  21(a-h).  Continued. 
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(a)  For  k  =  0. 
v  '  a 


(b)  Fork  =0.05. 
v  a 


Fig.  22(a-h).  Magnetic  moment  |  n^|  vs  T  for  k^  =  0.  70,  0  <  k^  <  0.50,  kQ  (as  indicated). 
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(c)  For  k  =0.10. 
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(d)  Fork  =0.15. 
a 


T(*K) 


Fig.  22(a-h).  Continued. 
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(f)  Fork  =0.25. 
'  a 


Fig.  22(a-h).  Continued. 
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(h)  For  kQ  =  0.35. 
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Fig.  22(p-h).  Continued. 
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(b)  Fork  =0.05. 
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Fig.  23(a-h).  Magnetic  moment  |  n^|  vs  T  for  k^  =  0.60,  0  ^  k^  ^  0.50,  k^  (as  indicated). 
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(c)  Fork  =0-10. 
x  a 
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(d)  Fork  =0.15. 
'  a 


Fig.  23(a-h).  Continued. 
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Fig.  23(a-h).  Continued. 
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(h)  For  kQ  =  0.35. 
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Fig.  23(a-h).  Continued. 
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(a)  For  k  =  0. 
'  9  a 


(b)  For  k  =  0. 05. 
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Fig.  24(a-h).  Magnetic  moment  |n^|  vs  T  for  kc  =  0. 50,  0  <  k^  <  0. 50,  kQ  (as  indicated). 
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(c)  Fork  =0.10. 
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(d)  For  kQ  =  0.15. 


Fig.  24(a-h).  Continued. 
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(e)  Fork  =0.20. 
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(0  Fork  =0.25. 
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Fig.  24(a-h).  Continued. 
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(g)  Fork  *  0. 30. 
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Fig.  24(a-h).  Continued. 
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(a)  For  kQ  =  0. 
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(b)  Fork  =0.05. 
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Fig.  25(a-h).  Magnetic  moment  (n^j  vs  T  for  k^  =  0.40,  0^  k ^-^0.50,  k^  (as  indicated). 
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(d)  For  kQ  =  0.15. 


Fig.  25(a-h).  Continued. 
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(e)  Fork  =  0.20. 
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(0  Fork  =0.25. 
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Fig.  25(a-h).  Continued. 
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(g)  Fork  =0.30. 
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Fig.  25(a-h).  Continued. 
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